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ABSTRACT 


Undesirable steady offsets result when a stationary, linear 
regulator using state feedback is subjected to constant disturbances 
and/or non- zero setpoints. To eliminate these offsets, the distur- 
bances and non- zero setpoints can be fed forward to the control. Only 
when the number of outputs is less than or equal to the number of 
control inputs can the outputs be maintained at arbitrary non- zero set- 
points. 

The constant disturbances may be estimated from the state using 
an observer. An alternative is to feed back integrals of the devia- 
tion in the outputs; this amounts to a special form of disturbance 
observer and has the advantage that the steady performance is insensi- 
tive to small deviations in the system model parameters. 

The state and the disturbance may be estimated using a constant 
gain Kalman filter. In this case we suggest that the constant distur- 
bances be modeled as exponentially correlated processes with long 
correlation times. As an alternative, a Kalman filter (neglecting 
disturbances) can be used to estimate the state and the estimated state 
and integrals of the measured output deviation are then fed back to 
the control. 

These results are applied to the problem of automatic landing of 
an aircraft in the presence of a steady crosswind. The goal of the 
control system is to maintain the lateral position and yaw attitude 
of the aircraft in alignment with the runway centerline. 
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.NOMENCLATURE 


For the context of general results: 


A Symmetric, positive semi-definite weighting matrix on the 

states for quadratic synthesis 

B Symmetric, positive definite weighting matrix on the controls 

C Matrix of control gains. Subscripts refer to associated 1 

vector, e.g, C , C , etc. 

x y 

F State dynamics matrix 

G Control distribution matrix 

H Measurement distribution matrix 

J Performance index for an optimization problem 

K Matrix of filter or observer gains 

L A special matrix useful in integral control synthesis 

M Matrix for feed-forward control 

P Steady estimation error covariance matrix 

Q Power spectral density of process noise 

R Power spectral density of measurement noise 

S Solution matrix of the Ricatti equation for quadratic synthesis 

sf Diagonal matrix of eigenvalues of Euler- Lag range system in 

quadratic synthesis. Subscript refers to sign of the real 
part. 

t Independent variable usually time or distance 

T Output distribution matrix 

u Control vector 

v Vector of integral control variables 

w Constant disturbance vector 

x State vector 

X Partition of the Euler- Lagrange eigenvector matrix corresponding 

to the state 

y Output vector 

z Measurement vector 

r Disturbance distribution matrix 

T) Impulse or white noise process disturbance vector 

X tfX Lagrange multiplier vectors 
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A 

u 

X 


(’) 
E( ) 

C) 

C) 

A 




Partition of the Euler- Lagrange eigenvector matrix corresponding 
to the adjoint variables 

Measurement noise vector 

Covariance matrix. Subscripts refer to associated vectors, 

e-g. X xu = e < xuT) 

Differentiation with respect to the independent variable 
Expectation operator (mean value) 

Estimate of a vector 

Error associated with the estimate of a vector, e.g. x=x-x 
Defined equal to 
Inverse of a matrix 

Transpose of a matrix (or column vector) 

Desired value of a vector 
Steady value of a vector 


For the context of specific results for the lateral motions of a DC-8 
aircraft in landing approach: 

V Reference airspeed of the aircraft (hft/sec , 1 hft = 100 ft) 

w Steady crosswind (ft/sec) 

x Independent variable = horizontal range (hft) 

y Lateral position (ft) 

6a Aileron (spoiler) deflection (deg) 

6r Rudder deflection (deg) 

€ Horizontal path azimuth angle (deg) 

cp Roll angle (deg) 

cp 1 Roll angle derivative (deg/hft) 

Yaw angle (deg) 

?jf T Yaw angle derivative (deg/hft) 

CO Turn rate (deg/sec) 
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Chapter 1 


Introduction 


The purpose of a regulator is to hold certain outputs of a system 
near desired set point values in the presence of disturbances* Two 
types of disturbance commonly encountered are (1) impulse (initial 
condition) disturbances and (2) constant (step) disturbances. To 
counteract impulse disturbances, feedback regulators are used; they 
are designed by time domain, frequency response, and state variable 
quadratic synthesis techniques (B-2). 

Constant disturbances arise from external influences on the system 
or from non-zero set points for the outputs. The effect. of constant 
disturbances on systems with regulators designed for impulse disturbances 
is to drive them through a transient to some steady offset. Compensation 
for this offset is often possible by feeding the set points and external 
constant disturbances forward to the control. Di Caprio and Wang, 1969 
(D-3), Belanger, 1970 (B-l), Kwakernaak and Si van, 1972 (K-2), and 
Power, 1973 (P-5) suggest breaking the non- zero set-point problem into 
two parts: (1) the determination of desired steady values of state and 

control which depend explicitly on the set points, and (2) the design of 
a regulator to control the deviations from the steady values. The similar 
problem of constant process disturbances is discussed in the scalar case 
by Athans, 1971 (A-2), and in the vector case by Johnson, 1971 (J-2). 
However, both of these treatments attempt to bring all of the state 
components to zero in the steady- state , which is usually not possible 
or desirable. In addition, the steady-state solution for zero output 
offset, when it exists, is often not unique. A static optimization 
procedure presented herein removes this ambiguity. 

If the constant disturbances are not measured directly, they may 
still be estimated (or otherwise accounted for) using only the available 
measurements . 

For single-input, single-output systems, a well-known compensation 
technique for constant disturbances is the addition of integral control 
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(C-l and D-l). This consists of feeding back the first integral of the 
deviation in the desired output. Through the application of frequency 
or time domain analysis, the integral control gain is adjusted for a 
suitable transient behavior. For a stable system, such integral feed- 
back causes the output being integrated to be driven to a zero steady 
offset . 

In the case of multi-input, multi-output systems, an observer (L-2 
and B-3) can often be designed to estimate the disturbance. The estimate 
is then fed to the control input. Johnson, 1971 (J-2) suggests this 
approach, although his formulation is somewhat different from the 
approach presented herein. 

The concept of integral control and its zero-offset properties can 
be extended to multi-variable systems. Johnson, 1968 (J-l), Latour, 

1971 (L-l), Athans, 1971 (A-2), and Anderson and Moore, 1971 (A-l) 
suggest modifying the performance index of the linear regulator, qua- 
dratic synthesis technique to include weighting on the derivative of the 
control input. After algebraic manipulation, the resulting control law 
is found to include integration of linear combinations of the state 
variables. One problem with this approach is that these combinations 
of state variables do not generally correspond to the desired outputs. 
Hence, the outputs will be driven to some offset unless the control can 
directly cancel the effect of the constant disturbance. This latter 
condition is overly restrictive. Porter, 1971 (P-1), Bryson, 1972 (B-5), 
and Kwakernaak and Si van, 1972 (K-2) suggest the more fruitful approach 
of adding to the performance index weights on the additional states 
formed by integration of the outputs. Meyer, Martin, and Power, 1973 
(M-l) present an interesting alternative formulation by weighting the 
output deviation and state derivative in the performance index. Some 
caution must be taken, however, since adding additional states may 
make a previously controllable system uncontrollable. Power and Porter, 
1970 (P-3) give the requisite controllability conditions for the case of 
additional integrated states. These conditions are found to be related 
to the ability to find some control that will yield zero offset. A 
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technique is presented herein which allows the use of the set-point 
gains for zero offset to augment a stable state feedback regulator with 
integrals of the output. The closed-loop poles of the overall system 
are those of the original regulator plus additional poles , corresponding 
to the integrations, with an arbitrary degree of stability. This 
technique requires very little additional computation and facilitates 
on-line switching from proportional to proportional-integral control 
modes. While preparing this report, reference was found to a similar 
pole assignment technique (P-4). The relation between the integral 
control and observer formulations is also explored. Davison and Smith, 
1971 (D-2) show that feedback of the state plus integrals of the output 
constitutes a minimal order realization of a system, operating on the 
state, for which the overall closed-loop poles can be assigned arbitarily 
and the steady output is zero. It is further shown here that integral 
control amounts to a special form of observer. 

In many cases, the full state is not known exactly. Instead, only 
certain noisy measurements are available. When an optimal filter (K-l, 
B-2) is used to estimate the state and constant disturbances, it is 
found that the asymptotic estimates of the constant disturbances are 
exact. Thus, the steady filter will have zero gains and cannot be used. 
To circumvent this problem, an observer design technique can be used. 
Kwatny, 1972 (K-3) found that by modeling the disturbance as a random 
walk process (integral of a white-noise process), such observers can 
effectively be designed. Alternatively, a filter which estimates the 
state alone (ignoring the constant disturbance) can be used in conjunc- 
tion with an integral control law, providing the desired outputs are 
directly measured. This latter technique has the advantage that the 
stationary output covariance does not depend upon the nature of the 
constant disturbance. 

Each of the above topics is discussed in detail in the remaining 
sections, and the results are applied to the problem of designing a 
lateral autopilot for automatic aircraft landing in the presence of a 
steady crosswind. 
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Chapter 2 


Review of Regulator Design for Impulse Disturbances 
with Exact Knowledge of the State 


2.1 - System Model 

The systems considered here can be represented by a set of constant 
coefficient linear differential equations 

x = Fx + Gu+IV) (2.1.1) 

where x = state vector 

u = control (input) vector 
7 ) = disturbance vector 
F = state dynamics matrix 
G = control distribution matrix 
r = disturbance distribution matrix 
(’) denotes derivative with respect to the 

independent variable (usually time or distance) 

The objective of regulator design is to keep the state, x, near 
zero in the presence of disturbances, 77 . Here we assume that (a) the 
disturbances are non-zero over periods of time short compared to 
response times of the system (impulse type disturbances), (b) the time 
average of the disturbances is zero (positive impulses are balanced by 
negative impulses), and (c) otherwise, the disturbances are unpredict- 
able. Since the disturbances are unpredictable except for the zero 
average, the best thing to do is to assume TJ =0 in the future and 
design the regulator to bring the state to zero from arbitrary initial 
conditions in an acceptably short period of time. Two techniques for 
designing such a control law are presented in the following sections. 

2.2 - Pole Assignment (B-3) 

For the linear system given by equation (2.1.1), let the control 


law be defined 



Thus, the closed-loop system will be given by 

x = (F-GC)x , x (t ) = x (2.2.2) 

o o 

The transient response will be determined by the eigenvalues 
(closed-loop poles) of the (F-GC) matrix. Stability is assured if all 
the eigenvalues have negative real parts. The dominant transient 
characteristics will be determined by the poles with the smallest 
magnitudes. If the system [F:G] is controllable, the eigenvalues can be 
arbitrarily assigned (complex values occurring in conjugate pairs) 
through the choice of C. To see this, consider the characteristic 
equation given by 

| si - F + GC| =0 (2.2.3) 

This determinant is a polynomial in the variable s and hence will have 
coefficients which are functions of the various elements of F, G, and 
C. Selection of the poles (roots of 2.2.3) prescribes these coefficients. 
Equating the coefficients yields algebraic (generally non-linear) equa- 
tions to be satisfied by the components of C. In the case of a single 
input system there are exactly enough components of C so the solution 
will be unique. However, for multi- input systems the solution for C 
is generally non-unique, leading to one of the difficulties with the 
method. 

Care must also be taken when assigning the magnitude of the eigen- 
values. Generally, if the magnitude of the closed-loop eigenvalues 
greatly exceeds the magnitude of the open-loop eigenvalues, large compo- 
nents of the control gain matrix (C) will be required. Large gains may 
lead to control saturation and noise amplification problems. 

2,3 - Quadratic Synthesis (B-2) 

Consider the following optimal regulator problem. We want to 
minimize with respect to control (u) the performance index 

CO 

J = 1/2 \ (x T Ax+u T Bu)dt (2.3.1) 

t 

o 
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subject to the differential equation constraint 

x = Fx + Gu - 

The matrices A and B are chosen subjectively to weigh the output devia- 
tion against control magnitude. The solution is the linear feedback law 

u = -Cx (2.3.2) 


where C = B 1 G T S , SF + F T S + A - SGB 1 G T S = 0 . (2.3.3) 

The solution of this non-linear Ricatti equation for the 

matrix (S) can be obtained efficiently by eigenvalue decomposition of 

the Euler-Language system (B~6). 


-1 

F , -GB G 


— 

X X 


i— 

X X 


O 



4- - 


+ — 


+ 

T 

-A , -F 


AA 


A A 


% 

1 

o 

*- — 


+ - 


+ - 

_ 


L + J 


(2.3.4) 


The matrix (S) is given in terms of the partitioned eigenvectors as 

S = Ax ” 1 (2.3.5) 

2,4 - Example: A Lateral Aircraft Autopilot 

Consider the design of an autopilot to align an aircraft with the 
runway in the final phase of the approach and landing. The lateral 
motions of an aircraft are well- approximated by a sixth order linear 
model (Appendices A&B). The model states are roll attitude (<p) and 
derivative (tp 1 ) yaw attitude (^f) and derivative (^ f ), horizontal path 
angle (€ ) , and lateral position (y) (See Fig. 1). The controls are the 
aileron (6a) and rudder (6r) deflections. The equations are written in 
terms of the horizontal range (x) as the independent variable, and the 
F,G matrices for a DC-8 aircraft in landing approach are given in 
Appendix B. 

Applying the quadratic synthesis technique described in the previous 
section, a performance index is chosen which weighs deviations in yaw 
attitude and lateral position against the control deflections. 
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Definition of Lateral Aircraft States 


VELOCITY w.r.L GROUNO 



HORIZONTAL 

sideslip w.r.t. the air (gives 
the aerodynamic forces) 

non-dimensional side wind (w/V) 

roll attitude 

yaw attitude 

flight path azimuth angle (y/V) 

lateral position 

range 

aileron-spoiler deflection 
rudder deflection 

FIG. 1 


P = 

Pw = 
9 = 

t|f = 

€ = 

y = 

X = 

6a = 
6r = 
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00 


J = 1/2 V [(y/a^ 2 + (f/a 2 ) 2 + ( 6a/b J 2 + (6r/b 2 ) 2 ]dx 

o (2.4.1) 

Choosing the weighting factors: 

a^ = 15 ft, a 2 = 2 de &» b 1 = 10 de S, and b 2 = 10 deg 

yields the feedback control law (see Appendix B): 

6a = -3.57143’ - 2.47 kp - 4.630f - 1.602?(r - 8.546C - 0.6592y -s 

1(2.4. 2) 

6r = -0.2328 Up’ - 0.4959p - 9.187^1" - 3.627^ - 2.917C - 0.0996yJ 

where: 6a, 6r, (p, ty, € are taken in units of deg. 
y is taken in units of ft. 

(p’ , f' are taken in units of deg/hft*. 

Using this control law, the closed-loop poles are shown in Fig. 2. 
The lateral position and yaw attitude are taken as zero when the air- 
craft is aligned with the extended runway centerline. The transient 
response to initial conditions of <p= 5 deg, $=5 deg, and y=15 ft 
is shown in Fig. 3. 


*lhft = 100 ft 
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DISPLACEMENT (ft) 


Response to Initial Condition Disturbance 
Using State Feedback Control 


LATERAL POSITION (y) 

20 

15 
10 



FIG. 3 
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Chapter 3 


Regulator Design for Constant Disturbances and Non-Zero Set Points 
with Exact Knowledge of the State 


- The Effect of Uncompensated Constant Disturbances 

Consider the linear system model of Chapter 2 with the addition of 
a constant disturbance vector and a specified output vector: 

x-Fx + Gu+Tw, y = Tx, (3.1.1) 

where w = constant disturbance vector, 

r = constant disturbance distribution matrix, 
y = output vector , 

T = output distribution matrix. 

The control objective is to bring and hold the outputs (y) to the set 

points (yj. 

a 

To illustrate the effect of a constant disturbance which is not 
compensated for, again consider the lateral aircraft control example. 
Using the feedback control law given before, the effect of a steady 
crosswind is shown in Fig. 4. The initial conditions are taken as zero, 
which means that the aircraft, initially in trim, suddenly encounters a 
crosswind which is then constant. The crosswind speed is taken as 25 ft/ 
sec 0=3 15 kts). Note that the lateral position (y), the yaw attitude 
(ffr), and roll attitude (<p) all go to appreciable steady offsets. 

3.2 - Design with Exact Knowledge of the Disturbances 

As a first step in the development, suppose that the disturbances 
(w) are known exactly. A reasonable control law would then be to feed- 
back the deviations in the state from a desired steady value (D-3, K-2, 
P-5). Thus, 

u = U g - C x (x-* g ) (3.2.1) 

The desired steady values u and x are chosen to satisfy the steady- 

s s 

state and zero-offset constraints and thus depend upon the set point and 
disturbances. 
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DISPLACEMENT (ft) 


Response to Step Crosswind Using State Feedback Only 



FIG. 4 




(3.2.2) 


Fx + Gu + Tw = 0 
s s 


Tx - y . .= 0 
s d 


(3.2.3) 


When the control law (3.2.1) is used in the system equations (3.1.1), 
the result is 


x = Fx + G[u -C (x-x )] + Fw 
S X s 


(3.2.4) 


or 


~ (x-x ) = F(x-x ) + G(u-u ) + Fx + Gu + Tw 
dt s s s s s 


= F(x~x ) + G(u-u ) 
s s 


(3.2.5) 


Thus, the problem of choosing the gain C is exactly that of the 

X 

regulator problem for Impulse disturbances of Chapter 2, where the 

deviations from steady values are used as state and control. The 

remaining problem is then the choice of appropriate steady values 

x and u 
s S 

The linear equations (3 . 2. 2) and (3 . 2. 3) have a solution if the rank 


of the matrix 
matrix 




is equal to the rank of the column- augmented 


F , G , Fw 

T > 0 ’ y d 


For arbitrary non-zero w and this requires that the number of 

outputs be less than or equal to the number of controls, and that 


the rank of 


F, G 
T, 0 


be equal to the number of rows. 


When the number of outputs equals the number of controls the 
solution is unique. 


X 


M 

, M 

s 



xw 

xy 

u 


M 

, M 

s 

wJ 


uw 

uy 


-r w 


(3.2.6) 
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where 



When the number of outputs is less than the number of controls, and 


the matrix 


r n 

F,G 

jr, o 


has full rank, an infinite number of solutions exist. 


A useful technique for specifying a solution is to minimize the performance 
index 


J = 1/2(x T Ax +u T Bu ) 
s s s s 


(3.2.8) 


subject to the steady-state and zero-offset constraints (3.2.2) and 
(3.2.3). This yields unique solutions of the same form as (3.2.6). 
Determination of the gain matrices in (3.2.6) for this case is discussed 
in Appendix C. 

Using (3.2.6) we may rewrite (3.2.1) in the form 

u = -c x + C y , - C w , (3.2.9) 

x yd w 9 


where 


C = M + C M 
y uy x xy 


C = (M +C M )r 
w uw x xw 


For the lateral aircraft control problem, there are two outputs, 
yaw attitude and lateral position, and two controls, rudder and aileron 
deflection. The choice of outputs was made among the three possibilities 
of lateral position, yaw attitude, and roll attitude. Obviously, the 
lateral position is of primary importance for an automatic landing. 
Deviations in roll atttitude will cause touchdown on one landing gear, 
while deviations in yaw attitude result in side skidding. The side 
skidding is deemed the more dangerous situation. 

Applying the previous results then yields the additional gains for 
non- zero set points and the steady crosswind. The control law becomes 
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(3.2.10) 


5a = -3.5710' - 2.4710 - 4.6301/r' - 1.602^ - 0.843td - 8.546€ 

- 0.6592(y-y d ) - 0.5732w 

6r = -0.23280* - 0.49590 - 9.187^' - 3.6221^ + 5.52td - 2.917€ 

- 0.0996(y-y ) + 0.4461W 

d 

where all quantities are as before except w which is the steady 
crosswind in units of ft/sec. 

The transient response to zero initial conditions is shown in 

Fig. 5 for y -ijr =0 . Note that both the lateral position and yaw 
a a 

attitude deviations go to zero in the steady state, while roll attitude 
maintains some offset. This offset is given as 

<p ~ (-0.108 deg-sec/ ft )w (3.2.11) 


3.3 - Design Using Estimates of the Disturbances 

The compensation method just described assumes perfect knowledge 
of the disturbances. If the disturbances cannot be measured directly 
it may still be possible to estimate them. One method of doing this is 
to synthesize a reduced order observer to estimate the disturbances 
(L-2, J-2). 


Since the state is assumed to be known exactly, its derivative, 
in principle, is also known exactly. Thus, let an estimate of w be 
defined as: 

w = k[ (F x+Gu+r*w) - x] , w(t ) - 0 (3,3.1) 


where K is a constant matrix to be chosen. The estimate error, 

_ A ^ 

w = w - w, is then given by 

w = kTw (3.3.2) 

If the matrix K is chosen so that kT is negative definite, the 
estimate error (w) will go to zero asymptotically, i.e. w -♦ w . 


The control used is the obvious modification of (3.2.9): 


u 


-C x + 
x 


C y - C w 
yd w 


(3.3.3) 
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DISPLACEMENT (ft) 


Controlled Response to Step Crosswind 
Using Exact Knowledge of State and Wind 



FIG. 5 


Since w = w+w , and w does not depend on x (Eq. 3,3.2) the 

eigenvalues of the controlled system are the eigenvalues of 

F-GC and of £ . 
x 

To avoid differentiating the state in (3.3.1), define 

A ^ 

= w + Kx ; then 
* ’ 



A 

w 


K(FxfGu+rw), 


w * " K* * 


(3.3.4) 


The deficiency of this method of compensation is that it depends 
upon precise knowledge of the system model. If the actual F, G, 
and r matrices differ from the estimated ones, a non-zero steady 
offset results. 


As an example, consider again the lateral aircraft control problem 

where the disturbance is the steady crosswind. Let the observer gains 

T —IT 

(K) be chosen so that K--CT (r T) _ r . Thus, the observer characteristic 
matrix (iT) becomes the scalar 

kT = -<T (3.3.5) 

and the observer pole is given by - a . In this case, (J was chosen 
as 0.3 hft ^ , which lies near the middle of the closed-loop regulator 
poles (see Fig. 2). The response to a 25 ft/sec crosswind is shown in 
Fig. 6. The steady values of lateral position (y) and yaw attitude (^) 
are still zero, but the transient is somewhat more pronounced than in 
the case of exact disturbance feedback. This is due to the lag in 
estimating the steady crosswind. 


3.4 - Design Using Estimates of the Output Deviations Produced 
by Uncompensated Constant Disturbances (Integral Control) 

An alternative to estimating the disturbances is to estimate the 
steady output deviations that would be produced by constant distur- 
bances if the system were not compensated. As we shall show, this leads 
directly to a generalization of classical integral control. 
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DISPLACEMENT (ft) 


Controlled Response to Step Crosswind 
Using Exact State and Estimated Wind 



FIG. 6 


Let 


u = -C x + C y 
x y o 


(3.4.1) 


Steady-state solutions, if any exist, must satisfy 

0 = (F-GC )x + GC y„ + Tw , 
x s y J o 

y . = Tx 
J d s 


(3.4.2) 


If the system is controllable, F-GC^ is negative-definite, so 


x = -(F-GC ) 1 (GC y +rw) 
s x y o 


(3.4.3) 


Substituting (3.4.3) into the second set of equations in (3,4.2) gives 

y d - iXw = LGC y y Q (3.4.4) 


where 


L = -T(F-GC )~ l . 
x 


Equation (3.4.4) is satisfied for all y, and w if 

d 


LGC - I , 

y 

y o = y d ■ lFw • 


(3.4.5) 


Gains can be found to satisfy (3.4.5) provided the number of outputs 

is less than or equal to the number of controls (and LG of full rank). 


IX -n 

An estimate of y = U. w can be obtained from 


w 


y = KL(x-Fx-Gu-]Tw ) 
w 

= -K(y -y ) 
w w 


(3.4.6) 

(3.4.7) 


Since y =0 , it follows from (3.4.7) that 
w 

y = -Ky where y 4 y - y 
w J w w J w 

Thus if K is chosen to be positive-definite y -+ 0 

w 

(3.4.1) into (3.4.6) yields 


(3.4.8) 

Substituting 


y = KLx - KL[ (F-GC )x + GCy - GC w + Tw] 
w x y o w 


(3.4.9) 
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(3.4.10) 


Using (3.4.5) and L(F-GC x > = -T , Eqn. (3.4.9) becomes 


A 



To avoid differentiating x 


KLx - K(-Tx+y d ) 

A A 

define y^ = y 
J * w 


KLx 


then 


where 


y. - w»-y d > ■ y»«„> = -n*« 0 > . ^ 

f. = f. * ■ ) 


u 


~C x 

X 


4 * 


C 

y 




j 


(3.4, 11 ) 


Now (3.4.11) is obviously a generalization of "integral control" from the 
single-input f single-output case to the multi-input, multi-output case. 

In steady-state, y + =0 =» Tx y^ , even if F, G, T were not estimated 
precisely. 


Since (3.4. S) does not depend on x , it follows that the closed- 
loop eigenvalues of the system using (3.4.11) are the eigenvalues of 
F-GC^ plus the eigenvalues of -K , which are determined separately. 

A slightly more familiar form of integral control is obtained from 
(3.4.11) by introducing a vector v such that 

y* £ Kv (3.4.12) 

Then (3.4.11) may be written as: 


v = Tx - y , v(t ) = -Lx(t ) , 
a o o 

u = -(C +C KL)x + C y. - C Kv , 
x y yd y 


(3.4.13) 


where LGC = I . 

y 

The eigenvalues of the closed-loop system are still the eigenvalues of 
F-GC x and of -K . The estimated steady output deviations that would 
be produced by constant disturbances if the system were not 
compensated are 

y w = K(v+Lx) . (3.4.14) 
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If integral feedback is added to a system already designed with 

state feedback without modifying the state feedback gains (note 

is modified to C +C KL in (3.4.13)), the transient behavior of the 
x y 

resulting closed-loop system is usually degraded (occasionally the 
system is even made unstable). Obviously, it is possible to design the 
state feedback gains simultaneously with the selection of the integral 
feedback gains so that the transient behavior is acceptable. This can 
be done, for example, by quadratic synthesis (see Appendix D) . 

A satisfactory state-plus-integral control system can always be 
found provided: 


(a) 

The system 

(F,G) is controllable. 


(b) 

The number of outputs is less than or 
number of control inputs. 

equal to the 

(c) 

The matrix 

pF,G~j has full rank, or, 

TO 

equivalently 


the matrix 
outputs. 

L 1 > U J 

LG has rank equal to the 

number of 


The major advantage of integral control is its relative insensi- 
tivity to errors in the system model. If the actual F, G, and T 
matrices differ from those used in the design, the transient behavior 
and the steady-state values of x and u are affected. However, 
unless the system is actually destabilized, the steady-output offset 
will still be zero. 


Again, consider the lateral aircraft control example. Since it is 
desired to drive both the outputs, lateral position (y) and yaw attitude 
(t), to zero offset, the feed forward gains will be given by (3.4.5) as: 


C 5= 

y 


(LG) 


.659 , -.843 

.100 , 5.52 


(3.4.15) 


Thus, for K 




-i 


and a = 0.3 hft " the additional integral 
control poles will be located at the same place as the observer pole of 
Section 3.3, The resulting feedback control law is 


21 



5a 


6r 


-5.7960* - 4.5670 - 7.703f 

+ 0. 659y - 0.1978 

a 

0.12030*- 0.67330 - 14.01^* 
+ 5,52 f + 0. lOOy 


- 2.768^ - 22.21€ - 2.547y - 0.843t J 

a 

^(y-y d )dx - 0.2502 ^f-^dx 

- 6.995^ - 5.658C - 0.2750y 

- .0298 ^<y-y d )dx - 1.657 

(3.4.16) 


The transient response to a 25 ft/sec steady crosswind from zero initial 

conditions with ^T,=y =0 is shown in Fig. 7. The steady offset for 
a a 

lateral position (y) and yaw attitude (^) is again zero, and the tran- 
sient response is almost identical to that of the first order reduced 
observer of Section 3.3 where the disturbance itself was estimated. 
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DISPLACEMENT (ft) 


Controlled Response to Step Crosswind 
Using Exact State and Integral Feedback on y and ty 



FIG. 7 





Chapter 4 


Review of Regulator Design for White Noise Disturbances Using 

Estimates of the State 


4,1 - System Model 

It is exceptional to be able to measure all of the state variables 
accurately enough that they may be considered to be known exactly, A 
more usual situation is that (a) only one or more functions of the state 
variables are measured, (b) these measurements contain random errors, 

(c) the process is driven by random disturbances. In this case the 
system no longer has a steady equilibrium state as assumed in Chapters 
2 and 3, However, it may have a "statistical steady state" with a 
stationary probability density for the state and control variables. 

If the process and the measurements can be modeled adequately as 
Gauss-Markov processes, the probability densities are completely speci- 
fied by mean values and covariance matrices of the state and control 
variables. 

We shall assume a stationary, linear model of the form 

x = Fx + Gu + D 7 (4.1.1) 

z = Hx + V (4.1.2) 


where z = measurement vector, 

H = measurement distribution matrix, 
77 = process disturbance vector, 

V = measurement error vector. 


We shall also assume that 17 and V are independent purely random 


processes (white noise) with 

E[77(t)] = 0 , ^ 

e[t? (t>| T (r )] = q 6 (t-r ) , 

E[>(t)] = 0 , \ 

E[y(t)y T (f )] - r 6(t^r) , 

E[rj(t)i/ T (r)] = 0 , j 


(4.1.3) 
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where 


Q = process noise spectral density, 

R = measurement noise spectral density. 


4.2 - Observer-Controller by Pole Assignment 

If the process and measurement noise are considered negligibly 
small (i.e. Q=0 , R=0 ), it is possible to "reconstruct 1 ’ the state 

using an "observer" (L-2), of the form 

x = Fx + Gu + K(z-Hx) (4.2.1) 


where x = estimate of x ♦ If the state x is observable with the 
measurements z , a gain matrix K can be found to assign the eigen- 
values (poles) of the estimate error system arbitrarily, where 


and 


k = (F-KH)x , 



x ss estimate error. 


(4.2.2) 


This estimate, x, can then be used with the controller feedback 
gains, C , discussed in Section 2.2: 

u = -c£ (4.2.3) 


The eigenvalues of the closed- loop observer-controller are the eigen- 
values of F-GC (the controller eigenvalues) plus the eigenvalues of 
F-KH (the observer poles). 


4.3 - Filter-Controller by Quadratic Synthesis 

A useful technique for the design of control systems with additive 
white noise in the measurements and the process is quadratic synthesis 
coupled with optimal filtering. The system is modeled as a vector 
Gauss-Markov process (Eqns. (4 . 1 . 1)- (4. 1. 3 ) , and the performance index 
is the expected value, given the measurement, of an integral quadratic 
penalty function: 


MinJ » C E { (x T Ax+u T Bu)| z} dt 
u y 


(4.3.1) 
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The minimizing solution has been shown (P-2) to be the optimal 
deterministic controller using the maximum likelihood estimates 
(Kalman filter, K-l) of the system states. Thus, u=-Cx and 

x = Fx+Gu+K( z-Hx) 

where C = B _1 G T S ; SF+F T S+A-SGB -1 G T S = 0 (4.3.2) 

and K = PH T R _ 1 ; FP+PF T -tTcT T -PH T R -1 HP = 0 


The steady covariance matrix of the state variables is given by 


where 


A ip ^ 

X = E ( xx ) - X + p 

(F-GC)X + X(F-GC) T + KRK T =0 


(4.3.3) 


P is the steady covariance matrix of the error in the estimates (x) 
which is uncorrelated with the estimate (x). 


The covariance matrix x is ttie stochastic analog of the deter- 
ministic steady response. The steady output covariance is a measure of 
how accurately the system is being controlled on the average. In fact, 
the performance index is a weighted covariance trade off between the 
state (or output) and the control. 


4.4 - Example: Lateral Aircraft Autopilot 

Consider again the lateral aircraft control problem. Suppose that 
instead of perfect knowledge of the state, three noisy measurements of 
roll attitude, yaw attitude, and lateral position are available. In 
addition, suppose the system is disturbed by gusty winds. The wind gust 
effects will be modeled as two independent Gaussian white noise 
processes: the first is due to lateral gusts, and the second is due to 
the lateral gradient of the vertical gusts. The power spectral densi- 
ties of the process and measurement noise and the RMS values corre- 
sponding to a 50 ft correlation length along the flight path are given 
in Table 4.1. 
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RMS 

Error 

0.872 

0.651 

0.366 

0.472 

0. 169 

3.09 


~ 

RMS 

Response 

1.47 

3.31 

0.500 





3.17 


The two poles near the origin correspond to the estimation of the 
lateral position by combining the attitude and position measurements. 
The relatively large time constant reflects the heavy filtering of the 
noise in the lateral position measurement. Also, the RMS roll response 
is considerably larger than the RMS yaw or path response. This is due 
to the use of banked turns as the primary controlling influence on the 
lateral position. 
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The RMS response in Table 4.2 can be compared to the RMS response 
shown in Table 4.3 where exact knowledge of the state is assumed. 


Table 4,3 



<pt 

deg/hft 

9 

deg 

f 

deg/hft 

t 

deg 

€ 

deg 

y 

ft 

6a 

deg 

6r 

deg 

RMS Response 

1.022 

2. 11 

0.309 

0.536 

0.407 

4.07 

5. 07 

2.92 
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Chapter 5 


Regulator Dealgn for Constant Disturbances and Non- Zero Set 
Points Using Estimates of the State 


5,1- The Effect of Uncompensated Constant Disturbances 


Consider the Gauss-Markov system model of Chapter 4 with the 
addition of a known constant disturbance vector, w: 


x = Fs + Gu+rw + Tv 
o 1 ' 

z = Hx + V , 


w = 0 


(5.1.1) 


If the filter-controller of Section 4.3 is used for the system, a 
statistical steady state will still be reached but it will have a non- 
zero mean value: 


e£x] - (F-KH) _1 T o w , "| 

E[x] - (F-GCT^GC (F-KH) _1 -l]r o W , \ (5.1.2) 

e[u] - -C{e[x3 + e[X]} . ^ 

Suppose that the aircraft with the filter-controller of Section 4.4 
is subjected to a 25 ft/sec steady crosswind. Since there is no direct 
compensation for a steady crosswind, it will go to steady offsets of 
-0.693 deg in roll, -4.38 deg in yaw, and 24.1 ft in lateral position. 
Such large offsets in yaw attitude and lateral position are unacceptable. 

5,2 - Design with Exact Knowledge of the Disturbances 

If the disturbances, w, are known exactly, this knowledge can be 
used in the filter for estimating the state and in the feed forward to 
the controls, i.e. 

u = -Cx + C y y d -C w w , (5.2.1) 

where 

x = Fx + Gu+Tw + K(z-l6) . (5.2.2) 

o 

The feedback gains C x and K are determined as in Chapters 2 and 4, 
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as though y d =0 , w=0 ; the feed forward gains C y and C w are 

determined as in Section 3.2. This system eliminates the mean value 
offset described in Section 5.1. 

5.3 - Design. Using Estimates of the Disturbances 

The constant disturbances are seldom known accurately so they must 
usually be estimated along with the state variables from the available 
measurements. However, the straightforward approach has a pitfall; by 
straightforward approach, we mean optimal estimation considering the 
disturbances as additional state variables with the random bias model, 

w = 0 , 

lCw(t©)] = 0 , E[w(t Q )w T (t o )] = X w » (5.3.1) 


i.e. the disturbances are constant, unpredictable, and vary from sample 
to sample. (In Section 5.1, the disturbances were assumed constant, 
known, and the same for all samples.) The controller with the augmented 
filter becomes 


ust-Cx+Cy.-Cw , 
x yd w 

x = Fx + Gu + Tw + K (z~Hx) , 
o 1 


(5.3.2) 


w = K 2 (z-Hx) 


where - 

0 = 


0 = 


P u hV 1 , K 2 = p^hV 1 , 

FP 12 +r <, P 22- P ll HVlHP 12 • 


(5.3.3) 


0 


T T -1 
P 12 H R HP 


12 ' 


However, (5.3.3) indicates that Kg=0 since both P^ 2 


and 


P =0 
22 


This occurs because the constant disturbances, w, can be estimated 
exactly by a time-varying Kalman filter. The K 2 =0 steady filter is 
useless for estimating the disturbances. Two alternatives to using a 
time-varying filter are (1) to use an exponentially-correlated model of 
the disturbances or (2) to use integral control. A random walk model of 
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the disturbances can also be used, but the steady covariance of the 
controlled state may be unbounded even though the filter is well 
defined. 


Exponentially-Correlated Model of the Disturbances 


A useful model of the disturbances is the exponentially-correlated 
model with a time constant, f , long compared to characteristic times of 
the system being controlled: 


* = - f w + TJ q 


e[t? ] = o , e[« (t>/V>] = Q 5(t-r) 


(5.3.4) 


where 


I Q o ■ X™ 


E[w(t)w T (t)] 


= given constant. 


This is essentially a slowly-changing bias model, which is more realistic 

than the ’’constant forever” bias model of (5.3.1). Here we must make a 

choice of T , whereas in a random walk model we must make a choice of 

Q ; the choice of Q is usually more nebulous than the choice of T • 
o o 

The steady filter that results from this model is again given by 
(5.3.2) and (5.3*3) except that the last two equations in (5.3.3) are 
replaced by 


0 = 
0 = 


FP 12 " T P 12 + r o P 22 _ P 11 H R “*12 

I ^22‘Xww 5 + P L HVlHP 12 



(5.3.5) 


and the second equation of (5.3.2) is replaced by 


/S 

W 


- ~ W + K(z-Hx) . 
T 2 


(5.3.6) 


With this filter, the system is predicted to reach a statistical 
steady state, and the usual techniques of predicting covariances may be 
used. Precisely speaking, the mean value offset of Section 5.1 will not 
be entirely removed since this filter will not estimate a truly constant 
disturbance exactly. However, for large T the mean value offset will 
be negligibly small. 


Use of the exponentially-correlated model to deBign a constant 

gain filter for estimating the lateral aircraft motions yields the 

filter gains shown in Appendix F where we assumed VT = 12500 ft and 
2 

v = (25 ft/sec) . The closed loop poles of the estimate error 

^WW 

equations are those of Fig. 8 plus an additional pole at -0.3 hft 
The steady ms response using this filter is shown in Table 5.1. In 
addition, the RMS response for w=0 and the steady response for a 
25 ft/sec steady crosswind are shown for comparison purposes with the 
preceding and succeeding sections. 

Table 5.1 



cp' 

<P 


' f 

€ 

y 

6a 

6r 


deg/hft 

deg 

deg/hft 

deg 

deg 

ft 

deg 

deg 

RMS 

1.595 

4.51 

0.635 

1.297 

0.713 

8.15 

10.12 

13.04 

RMS (w=0) 

i 

1.523 




0.647 

7.42 

6.87 

4.87 

Offset 

0.0 

-2.71 

0.0 



0.728 

-7.45 

12.18 


The offsets for ^ and y are due to the small error in estimating the 
constant disturbance and are insignificant when compared to the fluctuat- 
ing component in any given sample. 
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5.4 - Design Using Modified Integral Control 


Consider now the use of integral control with state estimation. 
Suppose a steady filter for the system has been designed assuming no 
constant disturbances: 

x = Fx + Gu + K^z-Hx) . (5.4.1) 


we then introduce integral control in the following form (cf. Section 
3.4): 


u = -Cx + Cy. - C (v+Lx) , (5.4.2) 

x yd v 


where 

v = Tx - y , + K(z-Hx) , (5.4.3) 

a i 

and C , C , C , L are chosen as in Section 3.4, while K is a 
x y v z 

constant matrix yet to be chosen. It is straightforward to show that, 

for any K , the eigenvalues of the closed-loop system are the eigen- 

values of F-GC plus the eigenvalues of F-K H plus the eigenvalues 
x 1 

of -LGC , so that C can be chosen separately rather than simulta- 
v v 

neously with and * 

K can be chosen to produce zero mean value offsets in the outputs, 
2 

For w=constant it follows from (5.4,1) and (5.4.3) that 


E (x) - (F-K 1 H)" 1 T > w , 

E(Tx-yJ - (K H-T)E(x) 
a £ 


where x=x-x=error in the estimate of x . Thus 


E(T*-y d ) - (K 2 H-T)(F-K 1 H) 


V w . 


(5.4.4) 


(5.4.5) 


If K can be chosen such that 

( T-K 0 H) (F-K H)r = 0 , 
z to 


(5.4.6) 


then the mean value offsets of the output will be zero. 

For the important particular case where the outputs are measured, 

i . e. 

z = Tx + V , (5.4.7) 
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so that T s H , it is obvious that (5,4,6) is satisfied if 

Kg = I = identity matrix . (5.4.8) 

The integrals of (5.4.3) become then simply 

v = z - y. . (5.4.9) 

d 

In general, the Kg matrix can be chosen to insure zero mean 
offset if the rank of the matrix [ T(F-KjH)~*T o ] is equal to the rank 
of the augmented matrix [ T(F-K H) -1 r : H(F-F^H) _1 T ] . Usually, this 
means that either the space spanned by T is included in that spanned 
by H or the dimension of the measurement is greater than or equal to 
the dimension of the steady disturbance. In the former case, it is 
possible to find K so that K H=T , and the outputs (y) are measured 

Z Zi 

directly with only white noise measurement error (y). In the more 
general second case, at least as many independent measurements as inde- 
pendant constant disturbances are required. In this second case, the 
stationary performance is sensitive to modeling errors in F and T ^ , 
whereas, in the first case when K 2 H=T * onl y sensitivity is to the 
measurement distribution H * 

Expressions for the steady covariance matrix of x are given in 
Appendix G. 

The above results are applied to the lateral aircraft control 
problem, using the undisturbed steady filter of Appendix E, the integral 
control law of Section 3.4, and the direct integration of the lateral 
position and yaw attitude measurements. The steady MS response to the 
noisy gusts and measurements and the steady offsets for the constant 
25 ft/sec crosswind are shown in Table 5.2. 
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Table 5.2 



<p’ 

<p 

t[r f 


€ 

MM 

6a 

fir 


deg/ hit 

deg 

deg/hft 

deg 

deg 

D 

deg 

deg 

RMS 

1.83 

4 * 15 

0*602 

1.02 

0.765 

o 

00 

CO 

9. 91 

5.21 

Offset 

0 

-2*71 

0 

0 

0 

0 

“7.64 

+12.5 


Notice that the lateral position RMS response is slightly smaller than 
using the state-disturbance estimator at the cost of increased RMS 
aileron control. 
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Chapter 6 


Conclusions 


Several techniques for constant disturbance compensation have been 
presented. The concept of integral control was generalized to the 
multi-input, multi-output output case. A procedure for synthesizing the 
integral control gains separately from the standard regulator-filter 
was presented. Using this technique, a regulator-filter whose transient 
response is satisfactory can be augmented with integral control for 
improved steady behavior. The state estimate filter is left unchanged 
and the state feedback gains are modified in a straightforward manner. 

The results were applied to the lateral aircraft control problem for 
automatic landing in the presence of a steady crosswind. It was found 
that state plus integral feedback provided an excellent technique for 
reducing the steady offset, without sacrificing good transient behavior. 
The use of a filter to estimate the state from noisy measurements coupled 
with the state plus integral feedback control laws resulted in acceptable 
RMS response and zero steady offset. 

The integral control technique of constant disturbance compensation 
was compared to compensation using disturbance estimators. It was found 
that integral control amounts to a special form of estimation. However, 
the RMS response may be larger for the integral control formulation 
because no use is made of the assumed knowledge of the disturbance dis- 
tribution as in the conventional estimator. This disadvantage of 
integral control must be weighed against the advantage that the steady 
offset remains zero even in the presence of modeling errors. The 
integral control technique does, however, exhibit sensitivity to bias 
or other long correlation disturbances in the measurements. Often, these 
measurement bias errors are not observable (as in the case of the lateral 
position measurement), and the conventional estimator does no better. 

In conclusion, then, integral control is likely to be most successful 
in situations where there are system model uncertainties and the measure- 
ments are unbiased. 
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Appendix A 


State Dynamics Matrix for Lateral Motions of a DC-8 Aircraft 
in Landing Approach with Open Loop Eigensystem 


ORDER OF SYSTEM • 6 


OPEN LOOP DYNAMICS MATRIX.... 


3* 9040-01 

0.0 

2.502D-G1 

2. 2400-01 

-2. 2400-01 

0.0 

<P r 

1.0000 00 

0.0 

0.0 

0.0 

0.0 

0.0 

<P 

5.0930-02 

0.0 

-1.087D-01 

-1.2760-01 

1.2760-01 

0.0 

r 

0.0 

0.0 

1,0000 00 

0.0 

0,0 

0.0 

t 

0,0 

5.4260-02 

0.0 

4. 5720-02 

-4.5720-02 

0.0 

c 

0.0 

0.0 

0.0 

0.0 

1.7450 00 

0,0 

y 


OPEN LOOP EIGENVALUES AND EIGENVECTORS.. (mf 1 ) 


COMPLEX EIGENVALUE! 1) 

< -0. 04483468) ( 0. 40648859) 


REAL EIGENVALUE ( 1) 

( 0.0 ).J! 0.0 ) 


REAL EIGENVALUE ( 2) 

( -0. 46048882). J! 0.0 ) 


COMPLEX EIGENVECTOR! 1) 

( -0.04483468)+JI 0.40648859) 

( 1.00000000.J! 0.0 ) 

( -0.262S0SS3).J! -0.20808545) 
( -0.43S38417)*J! 0.69380997) 

( 0.07822031)+J( -0.08434428) 

< -0,39431806 )>J! -0.29229680) 

REAL EIGENVECTOR! 1).. 

( 0.0 ) 

( 0.0 ) 

! 0.0 ) 

! 0.0 ) 

! 0.0 ) 

I 1.00000000) 

REAL EIGENVECTOR! 2) 

( -0.37782160) 

! 0.82047942) 

{ -0.00841106) 

{ 0.01826551) 

( -0.10934841) 

( 0.41437050) 


REAL EIGENVALUE ! 3) REAL EIGENVECTOR! 3) 

( 0. 00533819 ).J! 0.0 ) ( 0.00000167) 


! 0.00031376) 

! 0.00001625) 

! 0.00304391) 

! 0.00305910) 

! 0.99999064) 

REAL EIGENVALUE ! 4)...,.... REAL EIGENVECTOR! 4) 


( -0.00000000)+J! 0.0 


( 0 . 00000000 ) 
( 0 . 00000000 ) 
( 0 . 00000000 ) 
! - 0 , 00000000 ) 
! - 0 . 00000000 ) 
! 1 . 00000000 ) 
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Appendix B 


Control Distribution Matrix and State Feedback Gains 
for Lateral Motions of a DC-8 Airplane in Landing Approach, 
with Closed-Loop Eigensystem 


ORDER OF SYSTEM » 6 

NUMBER OF CONTROLS - 2 


OPEN LOOP DYNAMICS MATRIX.... 


-3* 904D-01 

0.0 

2. 5020*01 

2. 2400-01 

-2. 2400-01 

0.0 


1.0000 00 

0.0 

0.0 ' 

0.0 

0.0 

0.0 

CP 

*5.0930-02 

0.0 

*1.0870*01 

-1.2760-01 

1.276D-01 

0.0 

t 

0.0 

0.0 

1.0000 00 

0.0 

0,0 

0.0 

t 

0.0 

5,4260*02 

0.0 

4 . 5 720*02 

-4,5720-02 

0.0 

< 

0.0 

0.0 

0.0 

0.0 

1.7450 00 

0.0 

y 

STATE COST 

MATRIX 






0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0.0 

0,0 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

2.5000-01 

0.0 

0.0 


0 . 0 

0.0 

0.0 

0.0 

0.0 

0,0 


0.0 

0.0 

0,0 

0.0 

0.0 

4. 4440 

-03 


THE CONTROL DISTRIBUTION MATRIX.... 
Oa 5r 

1.2210*01 -3.057D-02 

0.0 0.0 

t.971D*03 6. S57D-02 

0.0 0.0 

0.0 -9. 772D-03 

0.0 0.0 


THE CONTROL WEIGHT I NO MATRIX 


1.0000*02 0.0 

0.0 1.0000*02 
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Appendix B (Cont. ) 


EIGENSYSTEM OF OPTIMAL CLOSED LOOP SYSTEM. . (lift’ 1 ) 


COMPLEX EIGENVALUE! 1).. ...... 

( -0.3230003D-J! 0.5024661B) 


COMPLEX EIGENVALUE! 2) 

( -0. 1199694U+J! 0.24002028 ) 


REAL EIGENVALUE ( 1) 

( -0.45287098>*J! 0.1) ) 


REAL EIGENVALUE ( 2) 

< *0. 25151727) *J( 0.0 ) 


THE CONTROL GAINS ARE: 

9 * v r 

-fi« 3.57180771 2.47071429 4.63055479 

-fir 0. 23302354 0.49602371 9.18713550 


COMPLEX El GEN VECTOR ( 1) 

( 0.53588176)*J! -0.19436612) 

( -0,75883202)+J( -0.57870318) 
( -0.32S0003D+J! 0.50246618) 

( 1.00000000)+J( 0.0 ) 

( 0.007725R8)+iJ( -0.021S65G5) 

< -0.06S19990)+J{ 0.01508147) 

COMPLEX EIGENVECTOR! 2) 

( 0. 201863 74 )-*J( -0.04316881) 

( -0.49691354)+J{ -0.59265542) 
( -0.01739570)+J( 0.00442558) 

( 0.04373718)+J{ 0.05061481) 

! -0.06875038 )*J( .0.13754744) 
( 1.00000000)*J( 0.0 ) 

REAL EIGENVECTOR! 1) 

! -0.37032118) 

( 0.81771894) 

! -0.00349654) 

! 0.00772084) 

( -0.11067944) 

< 0.42646943) 

REAL EIGENVECTOR! 2) 

! 0.12740354) 

( -0.50653995) 

! -0.00854898) 

( 0.03398965) 

! 0.12155221) 

< -0.84331625) 


♦ I 7 

1.60213745 8.54435546 0.$S91S493 

3.62741330 2.91641691 0.09957298 
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Appendix C 


Steady-State Solutions when the Number of Outputs is Less than 
the Number of Control Inputs 


The problem is to find x and u to minimize 

s s 

J = ^ (x^Ax + u^Bu ) , (C.l) 

2 s s s s 

subject to the steady-state constraints 

0 - Fx + Gu + Tw , (C. 2) 

s s 

y Tx (C . 3 ) 

J d s 

T 

Adjoining (C.2) and (C.3) to (C.l) by Lagrange multipliers X and 
T 

fi respectively, it follows directly that the unique minimizing values 

of x and u can be obtained by solving the linear equations 
s s 


-IT ™ 

F , -GB G , 0 


X 


r-rl 



s 


w 

. _T _T 

-A , -F , T 


X 


0 

T , 0 ,0 




_ y d^ 


where 

u = -B -1 G 1 X . (C. 5 ) 

s 

As an example, consider a steady banked turn at constant altitude 

for the DC-8 of Appendices A and B. Only the state variables , 

A 

and p = C - w are of interest in this case and it is clear that 

4>* = o , *r» = ^ , 

where V is forward velocity and (0 is the turn rate. The relevant 
constraint equations (with w=0 ) are therefore: 
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F 


G 


0 = 


r ^ 

0 f .2502, -.224^ 

0 ,-.1087, .1276 

.05426 , -1.000 , -04572 


s 

ilr 1 

s 



+ 


t ^ ' 

1224 , -.03057 " 

.008971 , .06557 

0 ,-.009772 



(C. 6) 


“ = [o, 1, <9 


f S 

s 


where 


is in units of hft/sec. 


(C. 7 ) 


For this case we take the weighting matrices of (C.l) as 

A = 0 , B = I . (C.8) 

Note there is only one output, ^ , while there are two control inputs 

( 6a and 6r ). With these data, the choices of $ ,ty’ ,6a ,6r 

s s s s s 

that minimize (C.l) while satisfying (C.6) and (C.7), obtained by 
solving (C. 4 ) and (C.5), are: 


s 


19.21 


~6a ** 
s 

1.000 

3|> 

6r 

0.97 22 

1 

_ s_ 


-0.3127 
-0. 1914 


tO 
V ’ 


(C.9) 


It is interesting to compare (C.9) with a "truly banked turn" 

where the aerodynamic side force is specified to be zero, i.e. 

-.045723 -.009772 6r =0 (there are now two outputs and two inputs), 
s s 

For this case it is easy to show from (3.2.7) that 


r~ — 

<P 


18. 43' 

s 




— 

1.000 

s 



/s_ 


-0.7335 


CO 

V 



-2.529 


to 
V ’ 


3.431 


(C.10) 
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For the zero sideslip turn 


‘ 9 
T s 

Notice that both the "truly banked" and zero sideslip turns require 
substantially larger control deflections than the control minimized 
turn. The steady side specific force, a, experienced by a passenger 
is given as 


18.77 


r- 


- -1 

CO 

V 

6a 


-1.576 

1.000 

s 

6r 


1.873 

0.000 


L sj 


L J 


a - g(18.43 6 C/V - cp ) , 

s 

where cp is in radians, CO in rad/sec, and V in hft sec This 
is equivalent to tilts from vertical of only 4.2$ and 1.8$ of the 
"truly banked" steady bank angle for the control minimized and zero 
sideslip turns, respectively. 
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Appendix D 


Quadratic Synthesis of Integral Control Gains 


An alternative method for determining state plus integral feedback 
control gains is to use quadratic synthesis (P-1, K-2, B-5). We tempo- 
rarily neglect the constant process disturbances and non- zero setpoints 
and augment the system with the additional states (v) given by v=Tx , 
Adding weighting on v in the performance index, yields an optimal 
regulator problem: 

00 

Min J = 1/2 \ (x T A n x + v T A 0 v + u T Bu)dt (D. 1) 

u * 

o 

Subject to x = Fx + Gu 

v — Tx 

where is arbitrary. When the controllability conditions are 

satisfied, the solution is 

u = -b“ 1 G T S 11 x - B -1 G T S V (D.3) 

11 1 ^ 

where + F T S n - S^GB^Vs.^ + A 1 + S^T + T T S^ 2 = 0 , 

(F-GB‘ 1 G T S 11 ) T S 12 + T T S 22 = 0 , (D. 4 ) 

and A 2 - S^GB'Vs^ = 0 


x(t • ) = x 
o o 


v(t ) = V 
o o 


(D. 2) 


To properly choose the initial conditions (v q ) , the performance index 

(D.l) will be minimized with respect to v (M-l). The performance 

o 

index is given in terms of the initial conditions as 



Thus, the minimizing is given by 


(D.5) 
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(D.6) 


v 

o 


- S 22 S 12 X 0 


When this control law is implemented with a non-zero setpoint the only 
modification required for a zero-offset steady state is to let 

v = Tx - y (D.7) 

However, the response to a setpoint change will be faster if feed- forward 
of the setpoint is also included. The resulting control law is thus, 

u=-Cx+Cyd-Cv (D.8) 

x y J v 
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Appendix E 


Filter Gains for the Estimation of Lateral Motions from Noisy 
Attitude and Position Measurements with Gusty Wind Disturbances 
and the Eigensystem of the Estimation Error 


ORDER OF SYSTEM - 6 

HUMBER OF CONTROLS - 2 

NUMBER OF OBSERVATIONS « 3 

NUMBER OF PROCESS NOISE SOURCES * 2 


OPEN LOOP DYNAMICS MATRIX..., 


-3.904D-01 

0.0 

2 • 502D-01 

2.2400-01 

-2,2400-01 

0,0 

1.000D 00 

0,0 

0.0 

0.0 

0.0 

0.0 

-5 • 093D-02 

0.0 

-1. 087D-0 1 

-1.27SD-01 

1. 2760-01 

0.0 

0.0 

0.0 

1.000D 00 

0.0 

0.0 

0.0 

0.0 

5 , 4260-02 

0.0 

4 .5720-02 

-4.5720-02 

0.0 

0.0 

0.0 

0.0 

0.0 

1. 7450 00 

0.0 

THE CONTROL 

DISTRIBUTION 

MATRIX.... 




0 * 

6 r 





1. 224D-01 

-3. 057D-02 





0.0 

O.fl 





t. 9710-03 

6.557D-02 





0.0 

0.0 





0.0 

-9 . 772D-03 





0.0 

4.0 






STATE DISTURBANCE DISTRIBUTION MATRIX.... 


2.2400-01 

0.0 

-1.276D-01 

0.0 

*.5720-02 

0.0 


-3.. 904D-01 

0.0 

-5.093D-02 

0.0 

0.0 

0.0 


POWER SPECTRAL DENSITY - STATE NOISE.... 


5.537D 00 0.0 

0.0 5.S37D 00 


MEASUREMENT SCALIN0 MATRIX 


0.0 

0.0 

1.000D 00 


2. 5000-01 0.0 o.n 

0.0 2. 500D-01 0.0 

0.0 0.0 1 . 000 D 02 


0.0 

0.0 

0.0 


1.000D 00 
0.0 
0.0 


0.0 

0.0 

0.0 


0,0 

1.0000 00 

0,0 


0.0 

0.0 

0.0 


POWER SPECTRAL DENSITY - MEASUREMENT NOISE.. 


48 


-V -G 



Appendix E (Cont.) 


EIGENSYSTEM OF ESTIMATE ERROR ERUATION (hft -1 ) 


COMPLEX E I GENVALUEt 1) 

( -1.03865067)*J< 1.01531494) 


COMPLEX El GENVALUE( 2) 

( -0. 5259 7023)*J( 0.61790966) 


COMPLEX EIGENVALUEt 3) 

( -0. 04887570 )*J( 0.04832768) 


COMPLEX EIGENVECTOR 1 )...... 

( 1.00000000)+d( 0.0 ) 

( 0.44759741)+J( -0.69319355) 

C -0.0 09 32262) +J( 0.06015641) 

( 0.08193031)*J( 0. 01090606 ) 

( 0.03536834)*J< -0.01701054) 

( -0.03751940)*J< -0.03386997) 

COMPLEX EIGENVECTORC 2) 

( 0. 28232492 )*J( 0.17562874) 

< 0.23021303)+J( 0.02854616) 

( 0.37325066)*J( 0.6186673S) 

( 1.00000000)*J( 0.0 ) 

< -0,13180058)+J( -0.16416696) 

( -0.331S2040)*J{ 0.19161289) 

COMPLEX EIGENVECTORt 3) 

{ -0.00305548)*J( -0.00348359) 
( -0.00208325)+J( -0.00215971) 
{ 0.004312S3)+J( 0.00491622) 

< 0.00628471)+J( 0.00555069) 

( 0.02577604)+JC 0.02682738) 

( 1.00000000)*d( 0.0 ) 


THE COVARIANCE OF THE ESTIMATION ERROR 


0, 76068438 

0*35913348 

-0.01907670 

0*04071315 

0.03108149 

0*35913348 

0*42368247 

-0.02072033 

0*00773122 

0.03185824 

-0*01907670 

-0.02072033 

0.13382472 

0.09957324 

-0.04144187 

0*04071315 

0*00773122 

0.09957324 

0*22297129 

-0.03374479 

0*03108149 

0. 03185824 

-0.04144187 

-0.03374479 

0.02865678 

-0.04427190 

-0.04427190 
0.00745114 
*0* 00535122 
-0.06528608 
0. 26659257 
9.55581974 

0. 00745114 

-0.00535122 

-0*06528608 

0.26659257 


FILTER STEADY STATE 

GAINS. 


1.43653394 

0*16285259 

-0.00044272 

1*69472986 

0.03092487 

0.00007451 

-0.08288130 

0*39829297 

-0.00005351 

0.03092487 

0*89188514 

-0.00065286 

0.12743295 

-0. 13497917 

0.00266593 

0*02980456 

-0.26114433 

0.09555820 

STATE RMS RESPONSE 

1.46850863 
<p 3.30690313 
r 0.49982089 
t 0.97850027 
C 0.65134758 
y 7.43165997 


CONTROL RMS RESPONSE 

0. 6.30470433 

Or 3.17164437 
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Appendix F 


Filter Gains for the Estimation of the Lateral Motions 
Including an Exponentially Correlated Wind Disturbance 
and the Eigensystem of the Estimation Error 


ORDER OF SYSTEM » 7 

NUMBER OF CONTROLS • 2 

NUMBER OF OBSERVATIONS • 3 

NUMBER OF PROCESS NOISE SOURCES ■ 3 


OPEN LOOP DYNAMICS MATRIX.*.. 


•3. 904D-01 

0.0 

2.502D-01 

2.2400*01 

-2. 2400-01 

0.0 

2.2400-01 

1.000D 00 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-5 . 0930-02 

0,0 

-1.087D-01 

-1.276D-01 

1.27BD-01 

0 . 0 

-1. 276D-01 

0 . 0 

0.0 

1. 000D 00 

0.0 

0.0 

0.0 

0.0 

0 . 0 

5. 4260*02 

0.0 

4.5 720*02 

-4 .5720-02 

0.0 

4. 5720-02 

0. 0 

0 . 0 

0.0 

0.0 

1.745D 00 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

-8. Q00D-03 


THE CONTROL DISTRIBUTION MATRIX 

Or 


1.224D-01 

0.0 

B.971D-03 

0.0 

0.0 

0.0 

0.0 


-3.05 7D-02 

0.0 

6.557D-02 

0.0 

-9, 772D-03 
0.0 
0.0 


STATE DISTURBANCE DISTRIBUTION MATRIX.... 


2. 2400-01 

-3. 9040-01 

0.0 

0.0 

0.0 

0.0 

-1.276D-01 

-5. 093D-02 

0.0 

0.0 

0.0 

0.0 

4.5720-02 

0.0 

0.0 

0.0 

0,0 

0.0 

0.0 

0.0 

l.OOOD 00 

POWER SPECTRAL DENSITY 

- STATE NOISE 

S.537D 00 

0,0 

0.0 

0.0 

5,5370 00 

0.0 


TD ^ ■£*. 



Appendix F (Cont. ) 


MEASUREMENT SCALING MATRIX..... 


0* 0 

1.0000 00 

0.0 0.0 

0.0 

0.0 

0.0 

0.0 

b.o 

0.0 

0.0 1.0000 

0.0 0.0 

00 0.0 
0.0 

0,0 

1,0000 00 

POWER SPECTRAL DENSITY - 

MEASUREMENT NOISE. 

ft 


2.5000-01 

0,0 

0.0 



6.0 

2. 5000-01 

0.0 



0.0 

0.0 

1.000D 02 



EIGENSYSTEM 

OF ESTIMATE 

ERROR EOUATION 

(6ft' 1 ) 



COMPLEX EIGENVALUE! 1) 

( -1,03699681)*J! 1.01796661) 


COMPLEX EIGENVALUE! 2) 

( -0. 52025652)»d ( 0.65*65307) 


COMPLEX EIGENVALUE! 3) 

{ -0.06»57769)*J! 0.06*67963) 


REAL EIGENVALUE ( 1) 

( -0, 300 232*6 ).J( 0.0 ) 


COMPLEX EIGENVECTORC 1) 

( 1.00000000)*J( 0.0 ) 

{ 0.66290256 )♦,)( -0.66328835) 

{ -0.0623876*)*«K 0.01656159) 

( 0.03737167)*J( 0.03625782) 

( 0.0S620321)+J( -0.025*9666) 

( 0. 01927629 )*J{ -0.06113573) 

( 0.369260S2)*J< -0.109060*2) 

COMPLEX EIGENVECTOR! 2) 

( -0.0353*373)+J< 0.12658826) 

( 0.01S60911)*JC 0.12575096) 

( -0.59215S69)+J{ 0.07566230) 

( -0.391S79S0HJC 0.56039075) 

( O.H5M6S7W! -0.06030693) 

< 0.1770*061)>J( -0.27307653) 

{ 1.00000000)+J( 0,0 ) 

COMPLEX EIGENVECTOR! 3) 

( 0*. 00003173 )+J( -0.00008351) 

( -0,00015*96)*J( -0.00006096) 
( 0.00001112 )+J! 0.00017166) 

( 0.00133603HJ! 0.00009723) 

( 0.02779562)*J( 0.027*6137) 

( 1.00000000)*J( 0.0 ) 

{ 0.02S76761)*J! 0.02717536) 

REAL EIGENVECTOR! 1).... 

( -0.11905073) 

( -0.07629656) 

( 0.15231280) 

< 0.1*072602) 

( -0.05695609) 

< -0.0900*263) 

! -0.95567316) 
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Appendix F (Cont. ) 


THE COVARIANCE OF THE ESTIMATION 


0. 79933060 

0.37985823 

0.37985823 

0.43549315 

-0.06643974 

-0.04773003 

0.00063180 

-0. 01646479 

0.04932665 

0.04233659 

0.00665403 

0.03868896 

0. 30049013 

0. 15956794 

0.00665403 

0.30049013 

0.03868896 

0.15956794 

-0.07677372 

-0.37077524 

-0.14197009 

-0.32021174 

0.29923830 

0.14371967 

9.78637978 

0.41803514 

0.41803514 

2.49140724 

ILTER STEADY STATE 

GAINS 

<P 

T 

* 

n 

1.51943293 

0.00252720 

1.74197261 

-0.06585916 

-0.19092014 

0.62444535 

-0.06585916 

1.11523408 

0.16934636 

-0.22463705 

0.15475583 

-0.56788036 

0.63827175 

-1.28084698 


ERROR 


-0. 0664397* 0.00063180 

-0.04773003 -0.01646479 

0. 10623749 0.15611134 

0.15611134 0.27880852 

-0.06562507 -0.05615926 

-0.07677372 -0.14197009 

-0.37077524 -0.32021174 


0.00006654 

0.00038689 

-0.00076774 

-0.00141970 

0.00299238 

0.09786380 

0.00418035 


0.04932665 

0.04233659 

-0.06562507 

-0.05615925 

0.03821684 

0.29923830 

0.14371967 
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Appendix G 


The Steady Covariance Using Modified Integral Control 

Using the equations given in Section 5.4, the estimates and 
estimate errors are found to satisfy the equations 


(F-GC + ) , -GC v , -K X H 
T , 0 , -K 2 H 


0 , 0 , CF-^H) 
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K_ 


-r 
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(G.l) 


Since y^ and w are considered to be deterministic, they affect only 
the mean value. The steady covariance then satisfies 

0 = (F-GC )*aa + Xaa ( F-GC ) T - GC yI - C T G T 
* A XX XX * v'xv A XV V 


0 - 0'-ec»>Xj v - oc^ vv + X jjT t * KjEK^ 

0 = T\a + T T + K RKj 
A xv A xv 2 2 

0 = (F-K H)\ — + X~~(F-K 1 H) T + K RkJ + r 

1 XX A xx 1 1111 


(G.2) 


and 

where 

and 


A xx 


■ Xvx=° 


and 


X— = p 

A xx 


0 = FP + PF T - PH T R 1 HP + 
T -1 

K 1 = PH R 


Thus, it is seen that the integral control law does not affect the 
steady performance of the filter, i.e. the estimate error is still P 
and the covariance between the estimates and the estimate errors 
remains zero. Only the steady covariance of the state estimates is 
changed by using the integral control law. 


X = + P 

A xx xx 


<uu - C *XJS C J * C v*vv C v + C .*iv C J 


r* 


T T 


53 



